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Two particle azimuthal correlation functions are presented for charged hadrons produced in 
Au + Au collisions at RHIC {^/s^ = 200 GeV). The correlation functions indicate sizeable 
asymmetries and anisotropics. The trend of the asymmetries is compatible with the presence 
of emission patterns associated with mini-jets. The magnitude and the trend of the differential 
anisotropics V2(j)t) and V2{Npart), provide important model constraints. 

1. Introduction 

Differential measurement of azimuthal (A0) correlation functions for charged hadrons 
is arguably one of the most important experimental probes for high energy-density nuclear 
matter created at the Brookhaven Relativistic Heavy Ion Collider (RHIC). The anisotropy 
as well as the asymmetry of such correlation functions are predicted to; (i) serve as a 
"barometer" for collective transverse flow and hence a probe for the equation of state 
(EOS) IHJ^, (ii) provide important constraints for the density and effective energy loss of 
partons and (iii) provide invaluable information on nuclear saturation scales 0. 

Full exploitation of these correlation functions require a good understanding of the mecha- 
nistic origin of two-particle azimuthal correlations at RHIC energies. Detailed systematic 
measurements of such correlations constitute important steps in this endeavor. 



2. Analysis and Results 

The colliding Au beams {y/s^ = 200 GeV) used in these measurements have been 
provided by the RHIC. Charged particles from mimum-bias triggered events were detected 
in the fully instrumented east and west central arms of PHENIX . Each of these arms 
subtends 90° in azimuth 0, and ±0.35 units of pseudo-rapidity rj. The axial magnetic field 
of PHENIX (0.5 T) allowed for the tracking of particles with p > 0.2 GeV/c {Sp/p ~ 1%) 
in the fiducial volume of both arms. Good track quality was ensured via matching and 
veto cuts to outer detectors in each arm. The Zero Degree Calorimeters (ZDC), were used 
in conjunction with the Beam-Beam Counters (BBC), to provide off-line selections of a 
wide range of centralities expressed as a fraction of the total interaction cross section. 

Two-particle azimuthal correlation functions were constructed via the ratio of two dis- 
tributions C{A(f)) = Ncor{A(f>) / ^mix{A4>) ^ where Ncor{A(p) is the observed A0 distri- 

*for the full PHENIX Collaboration author list and acknowledgments, see Appendix "Collaborations" of 
this volume. 
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bution for charged particle pairs selected from the same event, and Nmix{^(t)) is the A0 
distribution for particle pairs selected from mixed events. Mixed events were obtained 
by randomly selecting each member of a particle pair from different events having similar 
centrality and vertex position. Two correlation functions were obtained for each py-range 
of interest. In the first, charged hadron pairs were formed by selecting both particles 
from a reference range PTR^f^ which excludes the pt range of interest (i.e. a reference 
correlation |^). In the second, hadron pairs were formed by selecting one member from 
the Pt range of interest, and the other from the reference range PTn^f (i.e. an assorted-p^ 
correlation 

[1 .0 < pT < 2.5 (GeV/c)] [0.50 < pT < 1 .0 (GeV/c)] 




Figure 1. Assorted-pT correlation functions for charged hadrons of 2.5 < px < 3.0 
GeV/c. Correlation functions are shown for 1.0 < PTr^/ < 2.5 GeV/c (left panel) and for 
0.5 < PTjief < 1-0 GeV/c (right panel). The solid line represents a fit to the correlation func- 
tion, see text. 

Figure shows representative differential assorted-pr correlation functions for an event 
centrality cut of 20-40%. The hadron pairs were formed by selecting one member from the 
range 2.5 < px < 3.0 GeV/c and the other from a non-overlapping reference range PT^^f 
An additional requirement that particle pairs detected in the same PHENIX arm have 
the same charge (++ or --) was also imposed to reduce [but not exclude] non-flow corre- 
lations. The right panel of Fig. |l] indicates an essentially symmetric correlation function 
for 0.5 < PTjief < l-O GeV/c. By contrast the correlation function for 1.0 < Pt^^j, < 2.5 
GeV/c (left panel) shows an asymmetry at small A0. This asymmetry grows if narrow 
pseudo-rapidity cuts {Arj < 0.35) are imposed or if particle pairs are required to have dis- 
similar charge (+-). These trends are compatible with the effects of jet fragmentation [^]. 

The anisotropy and asymmetry of the generated correlation functions have been charac- 
terized via fits. A good representation of the data for the reference-pT and the assorted-pr 
correlation functions is obtained with the fit functions 

C(A(/)) = A ■ ea;p(-O.5(A0/(T)2) + ai(l + 2v^^cos{2A^)) 
C(A0) = A ■ exp{-0.5{A(j)/af) + ai(l + 2t;2, cos ( 2 A0)) 

respectively. Here the Gaussian term and the cos{2A(f)) term is used to characterize the 
asymmetry [at small Acp] and the anisotropy respectively. The anisotropy parameter for 
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the pt range of interest ^2, is given by v'2 = f2„/f2^ where f2„ and are the anisotropies 
extracted from the assorted and the reference correlations respectively. The value of 
can be compared to the second Fourier coefficient V2, commonly used to quantify the 
anisotropy with respect to the reaction plane |T0[. Differential values of V2 (hereafter 
referred to as V2) are summarized in Figs. |^ and |[ 
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Figure 2. V2 ys. pt for minimum bias events (left panel) and for several different centralities 
as indicated (right panel). Model comparisons are indicated for the minimum bias data. 
The inset shows model comparisons I^Q to data for the centrality cut of 20-40%. 

Figure ^ shows differential anisotropies V2{pt), which increase with pt followed by an 
apparent saturation for pt ^ 2.5 GeV/ c. This trend is similar for all centralities shown and 
is without doubt related to the collision dynamics. Hydrodynamic models provide 
a good description of the data for px up to ~ 1.5 GeV/c. However, these models do not 
predict the observed saturation of V2{pt) (see inset). The lack of saturation indicates that 
the model assumption that local equilibrium can be maintained until a sudden freeze-out 
hypersurface is reached, is invahd outside some finite domain of phase space. Results 
obtained from the jet dominated model HIJING, also show poor agreement with the data 
at high Pt- Models incorporating only strong jet quenching ^ or jet quenching acting 
in concert with hydrodynamic expansion ||^ do not provide a good description of the 
high-pT- data as well. By contrast, results from a recent implementation of a covariant 
transport theory |^ gives a good representation of the data if either extremely high initial 
gluon density or very large parton-parton scattering cross sections are employed in the 
calculations (see left panel of Fig. H). Results from a saturation model can also account 
for the observed anisotropies (see left panel and inset of Fig. ^ . 

The mechanistic origin of two-particle correlations in the transport and saturation mod- 
els are very different. Consequently, additional constraints are required to facilitate a 
distinction between them. One such constraint is the centrality dependence of V2- The 
saturation model predicts that V2 should scale as 1 / V Npart [§], where Npart represents 
the number of participants in the collision. For large opacities (high initial gluon density 
or very large parton-parton scattering cross sections) one expects the transport model to 
exhibit essentially Npart scaling. Fig. | show initial tests for 1 / a/ Npart scaling ( right 
panel), and Npart scaling (left panel) for charged hadrons of 2.5 < px < 4.0 GeV/c. 
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Figure 3. f2 vs. Npart (left panel), and \/Npart, (right panel) for charged hadrons of 
2.5 < pt < 4.0. Error bars indicate statistical errors only. 



The data indicates better overall agreement with Npart scaling for the centrality range 
presented but does not exclude Npart for relatively central collisions. 

3. Summary 

In summary, two-particle correlation measurements indicate sizeable anisotropies and 
asymmetries. The trend of the asymmetries are compatible with the asymmetric emis- 
sion pattern of mini-jets. The anisotropy {v2) increases with decreasing centrality and 
increases with px up to pt ~ 2.5 GeV/c. For higher pts the anisotropy saturates for 
each centrality cut. This important feature is not reproduced by hydrodynamic models 
but is well reproduced by the covariant transport theory of Molnar et al. [Q, and the 
saturation model of Kovchegov et al. More detailed measurements of the centrality 
dependence of w 2 may allow one to distinguish between the two very different mechanisms 
for two-particle correlations, implied by the latter models. 
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